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Abstract. An abstract approximation framework ami convergence theory for the identification of first 
and second order nonlinear distributed parameter systems developed previously by the authors and 
reported on in detail elsewhere are summarised and discussed. The theory is based upon results for 
systems whose dynamics can be described by monotone operators in Hilbert space and an abstract 
approximation theorem for the resulting nonlinear evolution system. The application of the theory 
together with numerical evidence demonstrating the feasibility of the general approach are discussed in 
the context of the identification of a first order rptasl-llnear parabolic model for one dimensional heat 
ronduclion/maxs transport and the Identification of a nonlinear dissipation mechanism (i.e. damp'"*) in 
a second order one dimensional wave equation. Computational and implementations! considerations, in 
particular, with regard to supercomputing, are addressed. 
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1. Introduction 

In this paper we report on the results of our efforts in the area of approximation for the identi- 
fication of nonlinear distributed parameter systems. The central focus of the present paper is our 
essential features of the underlying theory. The theoretical basis (i.e. convergence analysis) for the 
computational results to be presented below has been treated in detail in two of our earlier papers, 
[3] and [!>]. In [3] and [5] we developed nonlinear analogs of the general abstract approximation 
framework for the identification of linear systems given by Banks and Ito in [2]. Inverse problems for 
first order nonlinear evolution equations are handled in [3], while [5] is concerned with second order 
systems. We note that in another earlier paper, [4], we have developed an approximation framework 
for the estimation of parameters in nonantonomous nonlinear distributed systems. Although the 
general approach taken in [4] differs somewhat from the ideas in [3] and [5], and from those to be 
discussed here, the results presented there are certainly related, and at present, remain our only 
means of dealing with either linear or nonlinear nonautonomous systems. 

In the next section we briefly review some abstract functional analytic existence, uniqueness, 
regularity, and approximation results for nonlinear evolution equations in Banach spaces. In section 3 
we consider first order systems, define the class of inverse problems and evolution systems with which 
we shall be dealing and sketch the relevant approximation and convergence theory. We t hen consider 
an example involving a quasi-linear model for heat conduction and present the results (inclosing those 
aspects related to supercomputing) of our computational study and numerical investigations. In the 
fourth section we treat inverse problems for second order systems. In particular we consider the 
estimation or identification of nonlinear damping or dissipation mechanisms in distributed parameter 
models for mechanical systems. A brief fifth section contains some concluding remarks. 

2. Nonlinear Evolution Equations in Banach 
Spaces-Existence, Uniqueness, and Approximation Results 

We consider quasi-autonomous, in general nonlinear, initial value problems of the form 


(2.1) 

x(/) + Ax(i) 3 /(<), 

(2.2) 

H 

O 

1! 

H 


set in a Banach space X with norm ||x. In (2.1), (2.2) above we assume that T > 0, x° t Dom(A) = 
{x eX : Ax t Li(0,T\ X), and that for some wen the operator A + uil : X —* 2 X is m- 

accretive. In other words, that for some w e II we have (i) |xi — x 2 |x < |(1 + Aw)(xj — x 2 ) + A(yi — 
y2)U for every xi,x 2 e Dom(A), y x e Ax x ,yi t Ax 2 , and A > 0, and (ii) 77(7 + A(/l + w/)) = 
Uj-tDcmt(A) + w7))x = X for some A > 0. We note that A + wI in-accretive implies that for 

each A > 0 the resolvent of A + w/ at A, J(A; A + w7) : X — ♦ X, a single valued, everywhere defined, 
nonexpansive, nonlinear operator on X, can be defined by J( A; A + w7) = (7 + A (A -f w7))~'. 

A nonlinear evolution system on a subset H C X is a two parameter family of nonlinear operators, 
{U(t,s) : 0 < s < t < T), on ft satisfying U(t,s)<p t ft,U(s,s)tp = ip and U (t , s)U (s , r)ip = U(t,r)<p 
for every ip e SI and 0<r<s<<<T with the mapping (s,<) -♦ U(t,s)ip continuous from the 
triangle A = {(s,<) : 0 < J < I < 7) C R ! into X for each tp e ft. A strongly continuous function 
x : [0,T] — ► X is said to be a strong solution to the initial value problem (2.1), (2.2) if x(0) = x°, 
it is absolutely continuous on compact subintervals of (0,T), differentiable almost everywhere, and 
satisfies /(<) — x(t) t Ax(t) for almost every t t (0,7 ). 

It can be shown (see [7], [9], and [3]) that under the assumptions on A, f, and x° made above, 
a unique nonlinear evolution system {t/(f,s) : 0 < s < t < T], on Dom(A) with the following 
properties can be constructed. 

(i) |7/(<,s)ys — U(t,s)ij>\x < e u ’( ,_ *)|y5 - ip\x , for all <p, ip t Dom(A) and 0 <s<t<T; 
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t 

(ii) |£/(s + <,«)<£> — U(r + t,r)<p\x <2/ e a, (‘~ T )|/(r + s)— /(r + r)| A -dr, for all <p e Dom(A) and 

o 

all t > 0 such that s + t,r + t <T; 

(iii) If the initial value problem (2.1), (2.3) has a strong solution x, then x(t) = U(t,s)x(s), for 
0 <s<t<T; 

The strongly continuous function x given by x(<) = U(t,0)x°,t e [0,T], is referred to as the unique 
mild, generalized, or integral (see [6]) solution to (2.1), (2.2). It is immediately clear from (iii) above 
that when the initial value problem (2.1), (2.2) admits a strong solution, it and the mild solution 
coincide. 

The abstract approximation result which will play a fundamental role in the discussions to follow, 
is given in Theorem 2.1 below. The theorem we state here is similar in spirit to other related 
approximation results for nonlinear evolution equations - in particular, those that can be found in 
[8] and [10]. The proof of Theorem 2.1 is given in [3]. In the statement of Theorem 2.1 and elsewhere, 
for sets S n ,n — 0, 1,2, ... we use the notation lim S n D So to mean that given any s 0 t So , there 

n— >00 

exists s„ e S n , n = 1,2,... for which lim s n = so- 

n— ► oo 

Theorem 2.1. For each n e Z + = {1,2,3,...} let X n be a closed linear subspace of X, let f n t 
L\ (0,T;X n ), and let A n : X„ —>2 Xn be an operator on X„ with A„ +w/ m-accretive. Suppose that 
there exists a function g e Li(0,T) for which \f n (t)\x < g(t), and that lim Dom(A n ) D Dom(A). 

n — ► oo 

Suppose further that 

(i) lim /„(<) = /(<), a.e. te(0,T), 

n— *oo 

and that 

(ii) lim J(\;A„+wI)<p n = J(\;A+uI)<p for each tp e X whenever <p n e X n with lim <p n = <p, 

n— *oo n— ►oo 

for some A > 0. 

Then if {U„(t, s) : 0 < s < t < T) is the evolution system on Dom(A„) generated by A n and we 
have 

lim U n (t,s)<p n = U(t,s)>p 

n— *00 

uniformly in s and t for ( s,t ) e A, for each <p e Dom{A) and <p n e Dom(A„) with lim <p„ = <p. 

n— ♦ oo 

3. The Identification of First Order Systems with Dynamics Governed 
by Monotone Operators on Hilbert Space 
Let H be a real Hilbert space with inner product denoted by < •,• > and corresponding norm 
| • |, and let V be a real reflexive Banach space with norm || • ||. We shall assume that V is densely 
and continuously embedded in H. It follows therefore that V «— *■ H c — ► V* and that there exists a 
constant p > 0 for which ||y>||. < for all <p e H and \<p\ < m||^||, for all <p e V where || • ||, 
denotes the usual uniform operator norm on V* . We shall also use < •,• > to denote the natural 
extension of the H inner product to the duality pairing between V and V*. Let Q and Z be metric 
spaces and let Q be a fixed, nonempty, sequentially compact subset of Q. 

For each q e Q let 21(g) : V — *• V* be a single valued, everywhere defined, hemicontinuous (see 
[6]), in general nonlinear, operator from V into V * satisfying the following conditions. 

(A) (Continuity): For each <p t V, the map q — ► 21 (q)<p is continuous from Q C Q into V* . 

(B) (Equi-V-montonicity): There exist constants w e R and or > 0, both independent of q e Q, 
for which 

<%(q)ip-%.(q)ip,<p-i>> +w\ip-rl>\ 2 > a\\ip - rj>\\ 2 , 
for every <p, e V. 
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(C) (Equi-boundedness): There exists a constant /? > 0, independent of q £ Q, for which 


P(9)¥>l|. < /?( IMI + 1), 


for every ip £ V . 

For each q £ Q define the operator A(q) : Dom(A(q)) C H —* H to be the restriction of the 
operator 21(g) to the set Dorn(A(q)) = {<p £ V : 21 (q)ip t H}. It can be shown (see [3]) that A(q) 
is densely defined (i.e. that Dom(A(q)) — H) and that A(q) + uil is m-accretive. Let T > 0 and for 
each q e Q let /(• ;q) £ L x {0,T\ H) and let t i°(q) t H. We assume that the mapping q — ► u°(g) 
is continuous from Q C Q into II and that the mapping q — ► f(t;q) is continuous from Q C Q into 
H for almost every t £ (0,T). Also, for every z e Z, let u — ► ^(ujz) be a continuous map from 
C(0,T;II) into R + . We consider the abstract parameter identification problem given by: 


(ID) Given observations z e Z , determine parameters q e Q which minimize the functional 

<%) = <E >(«(?);*) 

where for each q e Q u(q) = u( ■ ;q) is the mild solution to the initial value problem in H 
given by 


(3.1) u(t) + A(q)u(t) = f(t ; q), 0 <t <T 

(3.2) «(0) = u°(g). 

Recalling the discussion in section 2 and our remarks above, it is clear that for each q £ Q , A{q) 
and /(• ; q) generate a nonlinear evolution system {U(t,s',q) : 0 < s < t < T} on H with the mild 
solution to the initial value problem (3.1), (3.2) given by 

u(<; q) = U (t, 0; q)u°(q), 0 < t < T. 

We develop an abstract Galerkin based approximation theory for problem (ID). For each n = 
1,2,... let H n be a finite dimensional subspace of H with H n C V for all n. We let P n : H — *■ H n 
denote the orthogonal projection of H onto H n with respect to the inner product < > and we 

make the standing assumption 

(D) lim | \P n ip — y>|| = 0, for each <p £ V. 

n— ► oo 

It is clear that assumption (D) together with the dense and continuous embedding of V in H yield 
that the P n tend strongly to the identity on H as well, as n — > oo. 

For each q e Q we define the Galerkin approximation ^(tf), n = 1,2, ..., to A(q) in the usual man- 
ner. That is, for <p n £ H n we set A n (q)<p n = rp n where ipn is the unique element in H„ (guaranteed to 
exist by the Riesz Representation Theorem applied to H n ) satisfying < %(q)<f n ,Xn >—< ’/’n.Xn >. 
Xn £ H n . We set /„(• ;q) = P n f( ■ ;g) £ L x {Q,T\H n ) and u°(g) = P n u°(q) £ H n . With these def- 
initions, it is not difficult to argue that for each n = 1,2,... and each q £ Q, A n {q) and /„(- ; q) 
generate a nonlinear evolution system, {U n (t,s;q) : 0 < s < t < T) on H n . Thus we consider the 
sequence of approximating identification problems given by: 

(/£)„) Determine parameters q n £ Q which minimize 

<M?) = $(««(?);*) 

where for each q £ Q,u„(q) = u n ( ■ ;q) is the mild solution to the initial value problem in 
H n given by 
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(3.3) 


u n (t) + A n (q)u n (t) = f(t,q), 0 <t<T 


(3.4) u n (0) = u°(q). 

The mild solution, u n (q), to (3.3), (3.4) is given by u n (t;q) = U n (t,0\q)u°(q), 0 < t < T. 

We may summarize the existence and convergence theory for solutions to problem (7 D n ) given 
in [3] as follows. Let {g n }^Li be a sequence in Q with lim q n = q 0 e Q. Assumption (D) and the 

n — oo 

continuity of the map q — » u° (q) imply lim u°(g n ) = w°(?o) in H. Similarly we have lim f n (t',<ln) — 

n— ► oo n— * oo 

f(t-,qo) in H for almost every t e (0 ,T) with | /„(• ;? n )| dominated by an L\ function which is inde- 
pendent of n. Assumption (D) also implies that lim Dom(A„(q„)) = lim H n D H = Dom(A(qo )) , 

n— oo n — ♦ oo 

and that for each A > 0, lim J(X]A n (q„) + uil)ip n = 7(A;A(?o) + u>I)<p in H whenever <p n e H n 

n — * oo 

with lim <p n = <p £ H. Thus, Theorem 2.1 yields 

n — -oo 

(3.5) lim u n (q n ) = Yim U n (-,Q;q n )u° n (q„) = U(-,0;q o )u°(q o ) = u(q 0 ), 

n — ►oo n— *oo 

in C(0,T;H). Similar arguments can be used to demonstrate that if {? m }m=i C Q with lim 9m — 9oi 
then 

(3.6) lim u n (q m ) = lim U n (-,0;q m )u°(q m ) = U„(-,0;q o )u°(q o ) = u n (q 0 ) 

m— *oo m-*oo 

in C(0,T;H) for each n = 1,2,... 

The compactness of Q, the continuity assumption on <$, and (3.6) are sufficient to conclude that 
for each n = 1,2,..., problem ( ID „ ) admits a solution q n e Q. Since C Q and Q compact, 

a convergent subsequence, {?n,}?Li> ma y be extracted from {9n}^i- If 9 = lim q nj , then for each 

3 J j—*OQ 3 

q £ Q (3.5) implies 


$(9) = 4H«(9 );z) = $( Jim u„ . (? n );z) 

J — CO 

= lim <Ku„,(9 ni );z) = .1™ ®n,(q nj ) 

j — * 00 * J } — oo 

(3.7) < .lim (9) = lim $(u„(g);z) 

J-+OO J — * OO 

= $( lim u n .(q);z) = $(u(g);z) 

j— 00 

= ®(9) 

and consequently that q is a solution to problem (ID). 

When the admissible parameter set Q is also infinite dimensional (when, for example, as is fre- 
quently the case, the unknown parameters to be identified are elements in a function space) it must 
be discretized as well. When this is in fact the case, the theory presented above requires the fol- 
lowing modification. For each m = 1,2,... let I™ : Q C Q — • ► Q be a continuous map with range 
Q™ _ l m (Q\ i n a finite dimensional space and with the property that lim I m (q ) = 9, uniformly on 

m — ► oo 

Q. We consider the doubly indexed sequence of approximating identification problems (ID™) where 
for each n and m (ID™) is defined to be problem (ID n ) with Q replaced by Q™. It can be shown 
that each of these problems admits a solution q ™ t Q™, and that the sequence {9^} will have a 
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Q-convergent subsequence, {q™f}, with limit in Q. If lim q^f = q e Q, then q can be shown to be 

j,k — oo 

a solution to problem (ID). The problems (ID™) involve the minimization of functionals over com- 
pact subsets of Euclidean space subject to finite dimensional state space constraints. Consequently 
they may be solved using standard computational techniques. 

We illustrate the application of our theory with an example involving the identification of a quasi- 
linear model for one dimensional heat flow (see [13], [14]). We consider the quasi-linear parabolic 
partial differential equation 


du 

dt 


(t,x) 



= f(t,x), 


t > 0,0 < x < 1, 


together with the Dirichlet boundary conditions 


u(f,0) =0, u(f,l) = 0, t >0 


and initial data 

u(0,x) = u°(x), 0 < x < 1. 

We assume that u° £ L 2 (0,1),/ £ Li(0, T; L 2 (0, 1)), and that q £ Cb(R), the space consisting of all 
bounded continuous functions defined on the entire real line and endowed with the usual supremum 
metric which we shall denote by doo(-, •)• We assume further that q satisfies 

(3.8) (q(0)0 - q(y)r])(0 - q) > a\0 - rj\ 2 , 0, t] £ R 

for some a > 0. (We note that if q is differentiable on R, then the Mean Value Theorem implies 
that the condition g'(0)0 + g(0) > oc > 0, 0 £ R, is sufficient to conclude that condition (3.8) holds.) 
To apply our framework we set H — L 2 (0,1) endowed with the standard inner product and norm, 

\D<j>(x)\ 2 dx)'j . In this case we have 

V* = H~ l ( 0, 1) and V ♦ H V* with the embeddings dense and continuous (n = 1). We take 
Q = Cb(R),^ = C(0,T; L 2 (0, 1)), and for given fixed values of «o , po , no , #o > 0 we take Q to be 
the Q-closure of the set 

{q £ Cb(R) : q(0) = q(~0), |?(0)| < Po, k'(0)| < co, q'(0)0 + q(O) > ao, for 0 t R, q(0) = 

constant for |0| > 0 q for some numbers 0 q satisfying O<0 q <0o}- 

A straight forward application of the Arzela-Ascoli Theorem reveals that Q is a sequentially compact 
subset of Cb(R)- If for each q e Q vie define the operator 21(g) : V —* V* by 

l 

<21 (q)<p,r/>>= j q(D<p(x))D<p(x)Dil>(x)dx, £ V, 
o 

then it is not difficult to show that assumptions (A), (B), and (C) are satisfied. Let {t,}f =1 with 
0 < <i < t 2 • • • < t„ < T be given, and for each z e Z define the least squares performance index 
#:C(0,T;L 2 (0,1))^R + by 


and set V — /fg(0,l) with norm || ■ || given by 


IWI = (/ 


(3.9) 



|t;(<i,x) - z(ti,x)\ 2 dx. 
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We consider the parameter estimation problem (ID) with Q,$, 21(g),/, and u° as defined above. 

For each n = 1,2,... let H n = span {<#,}?!/ where is the j-th linear B-spline on the interval 
[0,1] defined with respect to the uniform mesh {0, l/n,2/n, ..., 1}. That is 


(3.10) 


&(*) = { 


( 0 

nx - j + 1 
j + 1 — nx 

l 0 


o < X < ^ 

3 - 1 


<*<k 


j = 1,2, ■ • • , n - 1. Clearly H n C V = Hq(0, l),n = 1,2, ... Let P„ : H -* H n denote the orthogonal 
projection of Z, 2 (0, 1) onto H n with respect to the usual L 2 -inner product. Standard approximation 
results for interpolatory splines (see [15]) can be used to argue that assumption (D) is satisfied. 

We discretize the admissible parameter set as follows. For me Z + and q e Q set 


(/ m g)(0) = 9 e R 

j = 0 

where the 3p™(- ;0 g ), j = 0,l,2,...,m are the standard linear B-splines on the interval [O,0 g ] defined 
with respect to the uniform mesh {O,0 g /m, 20 q /m , ..., 0 q } and then extended to a continuous function 
on the entire positive real line via ipj'(9',9 q ) = ip™(9 q -,0 q ), 6 > 6 q . Using the Peano Kernel Theorem 
(see [12]) it can be argued that 

doo(I m q,q) = sup\I m q-q\ < “<r 0 , q e Q, 

l rn 

and consequently that lim I m q = q in CWR), uniformly in q for q e Q. 

m— ►oo 

For q m e Q m = I m (Q) the finite dimensional initial value problem (3.3), (3.4) takes the form 


(3.11) 


M n w n (t) + K„(w„(t); q m ) = F n (t), 0 < t < T 


(3.12) 


M n w n (0) = w° 


R n 1 is given 


where u>„(t) e R n_1 , M n is the (n — 1) x (n — 1) - Gram matrix whose (ij)-th entry is given 
by M‘ J =< <p' n , qP n >,F n (t ) and w° are the (n — l)-vectors whose i-th elements are given by 
F^(t) =< fn(t, ),ip n > and w°' =< u°,<p' n >, respectively, and K n (- \q m ) : R n-X - 
by 

ng m (nt; 1 )u 1 — ng m (n{t; 2 — v 1 })-^ 2 — v 1 }, t = 1 
nq m (n{v x — i> ,-1 }){i>‘ — u‘ -1 } 

—nq m (n{v i+l - td}){t> <+1 - n f }, i = 2, 3, ..., n - 2 
nq m (n{v n ~ 1 - v"- 2 })^- 1 - v n ~ 2 } 

+ng m (— nu n-1 )t) n-1 , i = n — 1 


K(v;q m ) 


for v e R n x . 

m 

If q m e Q m is given by q m (0) = ( Ij' x Pj l (\9\t9i m )> f° r 9 £ R, solving the identification problem 

3 = 0 

(ID™) involves the choosing of parameters (q™, q™, ...,g™, 9 q ™) T from a compact subset of R m + 2 
so as to minimize the functioned $(g m ) = $(u n (q m )-,z) where u„(q m ) is given by u„(t,q m ) = 
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n — 1 

J2 t £ [OiTl w *fh w n(' <Q m ) solution to the initial value problem (3.11), (3.12) in 

; = i 

R n_1 corresponding to q m . 

In order to actually test our scheme we let q"(8) = a(l — ,5e~ be ) for 9 £ R with a — . 9 and 6 = .5 
(note that 9Dq"(9) + q*(9) > .45 for 9 £ R and consequently that condition (3.8) is satisfied by q *) 
and set z(t, x) = 5e~‘(x — x 3 ), x s[0, 1], t > 0. Then setting 



x e [0, 1],/ > 0, and u°(x) = 5(x — x 3 ),x e [0,1], we used our scheme to attempt to estimate q* 
based upon the observations {z(.5j, ■)}}= 1 - 

All integrals that had to be computed numerically (i.e. some of the Ln inner products and the 
integral appearing in the definition of the least-squares performance index (3.9)) were computed 
using a composite two point Gauss-Legendre quadrature rule on [0,1]. For each n and m the IMSL 
implementation (routine ZXSSQ) of the Levenberg-Marquardt algorithm, an iterative steepest de- 
scent/Newton’s method hybrid, was used to solve the finite dimensional nonlinear least-squares 
minimization problem (ID™). For a given choice of the parameters q m , the initial value problem 
(3.11), (3.12) was solved in each iteration using the IMSL routine DGEAR with the stiff option 
operative. As an initial guess for q" we took q°(9) — 1, for 9 £ R with 9 q o = 4. 

All computations were carried out on a Cray X-MP/48 at the San Diego Supercomputer Center. 
Standard coding techniques which permit optimal vectorization were used whenever possible. These 
included the nesting of loops with the largest ranges the deepest, and the separation of vectorizable 
and non-vectorizable code into different loops. In general, we observed that in the absence of 
vectorization, the Cray was able to run our codes in approximately 1% of the time that it took an IBM 
3081. Vectorization on the Cray then yielded an additional speed-up factor of 17. Representative 
results that we obtained for various values of n and m are shown in Figures 3.1 - 3.3 below. 

The CPU times on the Cray for these runs ranged from about 3 seconds for n = 8, m = 1 to about 
180 seconds for n = 20, m = 4. The value of the performance index was reduced from <F„(g°) « 10 -2 
to $n(<7[T) ~ 10~ 4 . We solved the problem (ID™) unconstrained. That is we did not enforce the 
constraints in the definition of Q which render it and / m (Q) = Q m compact. Thus it was not 
surprising that, as we have seen before in the case of linear system identification, for each n, the 
inherent ill-posedness of the problem of identifying functional parameters began to cause difficulties 
as m was increased (see [1]). We were, to a certain degree, able to mitigate these instability effects 
with the introduction of Tikhonov regularization (see [11]). However, at least from a qualitative 
point of view, this is probably unnecessary since we seem to obtain reasonably good estimates with 
m relatively small. It is worth noting that we have also tested our approach on the much simpler 
problem of identifying constant parameters (for example, the estimation of the parameters a and b in 
the definition of q‘). In these tests it performed superbly with convergence to the true values of the 
parameters as n — ► oo immediately apparent. Finally, we also tested our scheme using discrete or 
sampled rather than distributed observations in the spatial variable although strictly speaking these 
examples can not be treated with our theory. With measurements taken at only one spatial point, 
x = .58 (i.e. with the observations {z(.5j, .58)}j 2 x ), the scheme’s performance remained essentially 
unchanged from that observed with distributed observations. We note that the existing theory for 
the case of linear dynamics (see [2]) can handle spatially discrete measurements. An extension of 
these results to the nonlinear case is currently being investigated but at present remains an open 
question. 

4. The Identification of Nonlinear Damping in Second Order Systems 

In this section we consider the identification of nonlinear dissipation mechanisms in abstract infi- 
nite dimensional second order elastic systems. In our treatment below we assume that the 
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Figure 3.3 Legend: q * , ff , q° — x — x — x . 

stiffness operator is linear. However we note that a similar approach can be used to identify a 
nonlinear stiffness operator in the presence of linear damping. We are currently looking into the 
extension of our theory to systems which involve both nonlinear stiffness and damping. A more 
detailed presentation along with proofs for the theoretical results we summarize below can be found 
in ^ 

Let the spaces H, V, V * , Q, and Z, and the set Q be as they were defined for abstract first order 
systems at the beginning of section 3. For each q £ Q let the operator 21 (g) e C(V, V*) satisfy the 
following conditions: 

(Al) (Symmetry) For all <p, £ V < ^{q)y>,rl> >=< <p,%{q)rl> >; 

(A2) (Continuity) For each <p e V the map q — *• 21 (q)<p is continuous from Q C Q into V*\ 

(A3) (Equi-V-Coercivity) There exist constants u e R and a > 0, both independent of q e Q for 
which < 21 (q)<p,<p > +w|y>| 2 > a||y?|| 2 , for all <p e V and q € Q; 

(A4) (Equi-Boundedness) The operators 21 (g) are uniformly bounded in q for q e Q. That is, there 
exists a constant j3 > 0, independent of q t Q, for which ||2l(g)j|, < 0\\<p\\, for all <p e V; 

Also, for each q e Q let the operator 23 (q) : £)om(18(g)) C V — * 2 V ‘ satisfy the foiling conditions: 

(Bl) (Domain Uniformity) Dom(*8(g)) = Dom(2 3) is independent of q for q e Q, and 0 1 Dom(23); 

(B2) (Continuity) For each g? £ Dom(2 3) and ip(q) e ®(g)g? the map q ip{q) is continuous from 
Q C Q into V* ; 

(B3) (Maximal Monotonicity) For all (^ 2 ,^ 2 ) t 23 ? = {(>p,rp) e Vx V * : >p e Dom(% %),ipe 

23(g)y>} we have < tpi - 1 P 2 , V?i _ <P 2 > > 0 with 23 ? not properly contained in any other 
subset of V x V m for which this monotonicity condition holds; 

(B4) (Equi-Boundedness) The operators 23(g) map V-bounded subsets of Dom{\ B) into subsets 
of V * which are uniformly bounded in q for g e Q. That is, if S is a V-bounded subset of 
£)om(23), the set 23(g)S is V*-bounded, uniformly in q for g e Q. 

Let T > 0 and for each g £ Q let u°(g) £ V, u*(g) e H, and /(• ;g) £ Li(Q,T;H). We assume 
that the mappings q — ► u°(q),q —* u^g), and g — ► /(<; g) are continuous from Q C Q into V,H, and 
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H respectively, for almost every t £ (0 , T). For every z £ Z let ( u,v ) — + <$(u,u;z) be a continuous 
mapping from C(0,T; V x H) into R + . The identification problem, which we shall again denote by 
(ID), takes the form: 

(ID) Given observations z e Z, determine parameters q e Q which minimize the functional 

$(?) = $(«(?)> ufa);*) 

where for each q £ Q,u(q) = u( ■ ;g) is the mild solution to the initial' value problem 

(4.1) il(t) + 23(g)tt(<) + 2l(g)u(<) 3 f(t\q), 0 < t < T 

(4.2) u(0) = u°(q), u(0) = u'(q). 

To make the notion of a mild solution to a second order initial value problem of the form (4.1), 

(4.2) precise, we rely on a reformulation as an equivalent first order system in a product space 
and then apply the abstract theory outlined in section 2. For each q t Q define the Hilbert space 
Xq — V x H with inner product (•, ■ ) q given by 

(4.3) ((ipi, V>i), (<P 2 , i>z))q =< 2t(g)vn, ^2 > +w<y?i,v? 2 > + <t/’l,V’2>- 

We denote the corresponding induced norm on X q by | • | ? . We note that our assumptions on the 
operators 21(g) guarantee that (4.3) indeed defines an inner product on V x H and that the Banach 
spaces {X ? , | • | ? ) are norm equivalent, uniformly in q for g e Q, to the Banach space X = V x H 
endowed with the standard product topology induced by the norm \(<p, ^) |x = (IMI 2 + 

For each q £ Q define the operator ^4(g) : Dom(A(q)) C X q — ► 2 Xq by 

MtiiV’iP) = {-ip, {%(<!)? + ®(g)^} n H), 

for (<p,ip) £ Dom(A(q)) = £ V x H : rfr £ V, {21(g)g5 + 23(g)V>} fl H / <j>}. It can be 

shown (see [5]) that there exists a 7 > 0, independent of g e Q, for which the operator A(g) + 7 1 is 
m-accretive on Dom(A(q)) C X q for each q £ Q- Also, for each q £ Q define F(- ; g) £ L 1 (0, T\ X q ) by 
^(^9) = (0,/(<;g)), for almost every t £ (0,T) and set x°(g) = (u°(g), ^(g)) e Xj. It follows that 
for every q £ Q, A(q) and F{- ; g) generate a nonlinear evolution system {U (t, s; g) : 0 < s < t < T} 
on Dom(A(q)) C Xq satisfying conditions (i) - (iii) given in section 2. Henceforth we shall assume 
that x°(g) eZ?om(A(g)) for each g £ Q, and by a mild solution to the second order initial value 
problem (4.1), (4.2) we shall mean the V-continuous function u(- ;g) given by the first component 
of the Xq (or X)-continuous function x( • ;g) = U(-,0;q)x°(q). We shall take u(- ; g) to be the 
H-continuous second component of x(- ;g). 

We note that if assumption (B4) is strengthened to the condition that the operators 23(g) map 
H-bounded subsets into ^‘-bounded subsets, uniformly in q for g £ Q, it can be argued (see [5]) 
that Dom(A(q)) — X q = X. In addition, since conditions (A3) and (A4) imply that Dom (21(g)) = 
{ip £ V : 2t(g)y> e H} is dense in V (see [16]) it is clear that the operators A(q) will also be densely 
defined when the set {<p £ V : 23(g)g? e H} is dense in H. In particular this will in fact be the case 
for all of the standard linear dissipation mechanisms - for example, air (23 ~ I), so called structural 
(23 ~ 2li), and Kelvin-Voigt viscoelastic (23 ~ 21) damping. 

With the existence and uniqueness of mild solutions on X q now demonstrated for each q £ Q, the 
q-uniform norm equivalence of X q and X will allow us to subsequently ignore the q-dependence of the 
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state spaces and to develop our approximation theory and convergence results on the q-independent 
space X. 

Once'again, as was the case with first order systems, our approximation theory is of Galerkin 
type. For each n = 1,2,... let H„ be a finite dimensional subspace of II with H n C V. Let 
P n \ H —* H n denote the orthogonal projection of H onto H n with respect to the standard inner 
product on H, < >, and we assume that P n (Dom(*B)) C Dom(23), for all n. We also again 

assume that condition (D) given in section 3 is satisfied. For each n = 1,2,... and each q £ Q 
let 21 n (g) e C(H n ) and 23„(g) : Dom(23 n ) C — ► 2 Hn denote the Galerkin approximations to 

21 (g) and 23(g) respectively. That is for <p n e H n ,^- n (q) l Pn — 4'n where xp„ is the unique ele- 
ment in H n which satisfies < 21 (g)y>„,Xn >=< f n ,Xn >,Xn £ H n , and for ip„ e Dom( 23 n ) = 
Dom(tB) D H n , <3 n (q)ip n = {xp n e H n .< rp, Xn > = < ^n,Xn >,Xn £ H n for some ip £ 23(g) n }. 
We set /„(• \q) = P„f ( • \q) £ Li(0,T; H n ), u°(g) = P„u°(g), and u*(g) = P n u l (q), and consider 
the sequence of approximating parameter identification problems given by: 

( I D n ) Determine parameters q n £ Q which minimize the functional 

$n(<?) = $(u n (q),ii n (q)-,z) 

where for each q £ Q,u n (q) = u n (- ;q) is the mild solution to the initial value problem in H„ given 

by 

(4.4) u n (t) + <B n (q)u n (t) + 2 n (q)u n (() 3 f n (t;q), 0<t<T 

(4.5) u„(0) = u°(g), u„(0) = u\(q). 

We again use the theory in section 2 to define what is meant by a mild solution to the second 
order initial value problem (4.4), (4.5) in H„. For each n = 1,2,... let X n = Hn x //„, and for each 
q £ Q define the operator A n (q) : Dom(A n ) C X„ — *• 2 X ' by 

A n (q)(<p n ,xp n ) = (-^n, a n(9)y>r. + 

for (tpn.ipn) £ Dom(A n ) = H n x Dom(\ B n ). Set F „( • \q) = (0 , /„(• ;?)) £ Li(0,T]X n ) and 
x°(g) = (u°(g), u^(g)). We assume u\(q) £ Dom(*B„) so that x°(g) £Dom(A n ). We define the mild 
solution to (4.4), (4.5) to be the first component of the function x„( ■ ; q) = U n ( •, 0; q)x„(q) £ C(0,T; 
X n ) where {?/„(<, s; 9 ) : 0 < s < t < T) is the nonlinear evolution system on Dom(A n ) generated 
by A n (q) and F„(- ;g). That such an evolution system in fact exists can be argued as it was for the 
corresponding infinite dimensional second order system using the definitions of the operators A„(q) 
and 23 n (g), and the function /„(• ;q) (see [5]). The function ti n (- ;g) is obtained from the second 
component of x n (- ; q). 

By using condition (D) to argue resolvent convergence, i.e., that for each A > 0 sufficiently large, 
lim J(\;A n (q n )+u>I)(<p n , 1 Pn) = J(^',A(qo)+vI)(<j>,ip) in X for some u> £ R whenever (q>,ip) £ X 

n — ► 00 

and (<p n ,tp n ) £ X n with lim (<p n ,xp n ) = (y>, i>) and q n ,qo £ Q with lim q n = q 0 , we are able to 

n — 00 n— • 00 

apply Theorem 2.1 to obtain that lim u n (q„) = u(qo) in C(0,T; V) and lim u n (q n ) — w(<7o) in 

n — ► 00 n — ► 00 

C(0,T;H) whenever q n , qo £ Q with lim q„ = go- A continuous dependence result analogous to (3.6) 

can also be obtained in this fashion. Then using estimates in the spirit of those given in (3.7) we find 

that solutions q n to the problem ( ID „ ) exist and that the sequence {gn}^Li admits a Q-convergent 

subsequence, {gny}fLn with lim g n = g and g a solution to problem (ID). The discretization of 
3 3 n— *00 3 
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the admissible parameter set Q can be carried out, and a subsequent convergence theory established 
exactly as they were in section 3 for first order systems. A complete and detailed discussion of the 
results for second order inverse problems which we have summarized above can be found in [5]. 

To illustrate the application of our approach we consider an inverse problem involving the iden- 
tification of a nonlinear damping functional in a one dimensional wave equation. Let Q = Cs(R) 
endowed with the usual supremum norm, and let Q be the Q-closure of the set 

{q e Q - q(0) = -q(-9), 9q{9) > 0, for 9 e R, |g(0)| = q(9 0 ), 
for \8\>Q 0 ,q e H l (-9 0 ,9 0 ), |g|//i(-0 o ,0 o ) < Ko, 
q'(8) > 0, for a.e. 6 e (-9 0 ,9 0 )} 

where 9o and Ko are given positive constants. It is not difficult to show that Q is a compact subset 
of O. 

For each q £ Q we consider the one dimensional wave equation with nonlinear damping given by 

+ =/(<,*), <>0, 0 < x < 1 

with boundary and initial conditions 

u(<,0) = 0, ti(t,l) = 0, t > 0, 
u(0, x) = u°(x), ^ (0, x) = u ! (x), 0<x<l 

where E e Loo(0, 1) with E(x) > E 0 > 0, a.e. x £ (0,1),/ £ L 2 ((0,T) x (0,1)), u° £ Hq( 0,1), and 
it 1 £ L 2 (0, 1) are given. We set H = L 2 (G, 1), F = /fo(0, 1) and V * = Lf _1 (0,l). The operator 
21 £ C(V, V*) given by 

l 

<21 <p,i>>= J E(x)D<f{x)Dip(x)dx, <p,4> £ 
o 


is easily shown to satisfy conditions (Al) - (A4). For each q £ Q we define the operator 23 (q) : V — ► 
V* via 

l 

<23 (q)<p,if>= J q(<p{e))rf{9)d0, <p,4> £ 7fo(0, 1). 

o 

(Note that in this case we in fact have 72.(23 (q)) C H.) With the set Q as it has been defined above, 
it is clear that conditions (B1)-(B4) are satisfied and moreover that 23(g) maps H-bounded subsets 
of V into ^'-bounded subsets, uniformly in q for q £ Q. 

We take the observation space Z to be xJTj (R* x i 2 (0, 1)} and a weighted least-squares perfor- 
mance index, $, of the form 


(4.6) 


$(vmM) = 

*=1 


f 

- Z lj\ 2 + ffi / \HU,x) - zf{x)\‘ 
{ o 


dx 


for <p £ C(0,T\V),ip £ C(0,T]H), and z = ((z}, 2 ?), (z 2 ,z 2 ),. ..,(z^, 2 ^)) £ Z with p iy <Ti > 0, i = 
1,2, 0 < ti < t 2 < ■■ ■ < t v < T, and 0 < x x < x 2 < • • • < x t < 1. 
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As in our first order example, we employ linear spline based state approximation. For each n = 
1,2,... let H n = span{y>J l }"r 1 1 where the are the standard linear B-splines on the interval [0,1] 
defined with respect to the uniform mesh {0, 1/n, 2/n, ..., 1} as given by (3.10). Let P„ : H — ► H n 
denote the orthogonal projection of 1) onto H n with respect to the standard Z^-inner product. 
We again use linear interpolating splines to discretize the admissible parameter set Q. For m £ Z + 
and q £ Q set 


(4.7) ( i m g)(e ) = X>(jV™Wr( VW'gntf) 

i= i 

for 6 £ R, where the ^>T(-\9 q ) are as they were defined in section 3, and 9 q is that number in (0,#o) 
for which |</(0)| = q(9 q ), |#| > 9 q ■ (Note that in this case the lower limit of the sum in (4.7) is 1 
rather than 0 since q £ Q implies q(0) = 0.) We again have that condition (D) is satisfied and that 
lim I m q — q, uniformly in q for q £ Q. We set Q m = I m (Q). 

m— *00 

For H n and Q m as defined above, the finite dimensional initial value problem (4.4), (4.5) takes 
the form 

(4.8) M n w n (t) + C n (w n (t);q m ) + I< n w n (t) = F n (t), 0 < t < T, 

(4.9) M n w n (0) = w°, M n ii>„(0) = w £, 

where the (n — 1) x (n — 1) matrix M„ and the (n-l)-vectors f n (t) and u>° are as they were defined 
in section 3, tu* is the (n-l)-vector whose i-th component is given by w = < v},ip x n >, K n is 
the (n — 1) x (n — 1) matrix whose (ij)-th entry is given by K ^ =< Eip J,,^ >, and C„(- ;q m ) : 
R "- 1 R n_1 is given by 

J. 

A 

C£(v;g m ) = J {nx — 1 + 1} q m ({nr - :} { v 1 - v 1 *" 1 } + v‘) dx 

»-i 

n 


i± 1 

n 

+ / {. -4-1- nx} q m ({nx - *'} { v ,+1 - v' } 4- v') dx, 

jL 

n 

i — 1, 2, • • • , n — 1, for v e R n_1 with v°,v n = 0. If w„(- ; q m ) is the solution to the second order 

n — 1 

initial value problem (4.8), (4.9) corresponding to q m £ Q m , then u n (t;q m ) = w n(^>1 m ) l Pn an d 

;'=i 

«n(<;9 m )= E for t £[0,T]. If q m £ Q m is given by q m (9) = E 9jVj"(|0|; d q m ) sgn(0), 

;=i J=i 

9 eR, the identification problem (ID™) becomes one of determining parameters (gY 1 , ...,q%,9 q m) in 
some compact subset of R m+1 which minimize $ n (q m ) = $„(u„(q m ), u„(q m ); z). 

In order to actually test our scheme, we set 


f 0\9\ a sqn(e) -9 q . <0< 9 q . 

I P\8q'\ a sqn(6) |0| > 0 q . 
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with /? = .15, a = 2, and 0 q • = 2.5. With 


y(t,x) = {3cos(l-7rt) + 2sin(\-irt))sinnx 
o L 


for t > 0 and x £ [0, 1], we set 




with E = 1 , u°(x) = y(0,x) = Ssimrx, and ^(x) = |^( 0 , x) = 7rsin7rr, for t > 0 and x e [0,1]. 
For observations upon which to base our fit, we took z = {(z/j, 2 ?)}i£i with z}^ = y(.5i, .12) and 
zf(x) = a;), x e [ 0 , 1 ], * = 1 , 2 ,..., 10 . As an initial guess we set 


f .60 —1-2 < (9 < 1.2 

1 .6 sgn(6) |0| > 1.2. 


The weights {pi},-£i and in the performance index (4.6) were all set equal to 1. 

Using the same computational techniques and resources (both hardware and software) that we 
used for the first order example described in the previous section, we obtained the results plotted 
in Figures 4.1 - 4.3 below. The CPU times on the Cray X-MP/48 for these runs ranged from 84.96 
seconds for n = 8, m = 3 to 1032.78 seconds for n = 20, m = 3. When n = 20, the value of the 
performance index was reduced from (?°) ~ 6.0 x 10 -2 to ^(f^) « 2.0 x 10 -3 . For other 
values of n, the reduction in <£„ was less pronounced. In this particular example we found (and it is 
apparent from the figures) that truly satisfactory results could not be obtained until n was chosen 
sufficiently large. However, as is clear from Figure 4.3 the scheme performed extremely well when 
n = 20. Once again, as in the case of a first order system, we found that although our theory does 
not apply, similar results could be obtained using a performance index involving spatially discrete 
measurements of velocity. As expected, since the problems (ID™) were solved unconstrained (i.e. 
the compactness assumption on Q, and therefore Q m for each m, were not enforced) the presence of 
instabilities became apparent for each n with m sufficiently large. 


5. Concluding Remarks 

In this paper we have summarized the theoretical framework for the identification of nonlinear 
distributed parameter systems which we have developed elsewhere ([3] and [5]), and, more impor- 
tantly, have for the first time provided numerical evidence that our approach is indeed feasible and 
in fact performs well. In the case of second order systems, while our focus here has been on the 
identification of nonlinear damping in systems with linear stiffness, our theory is easily modified to 
handle the estimation of a nonlinear stiffness operator in the presence of linear damping. We are 
currently studying the extension of our results to second order systems which simultaneously involve 
both nonlinear stiffness and damping. Further numerical studies involving supercomputing are also 
presently underway. In addition to continuing our efforts using simulation data, we intend to test 
our schemes using experimental data in the near future. 
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identification of a first order quasi-linear parabolic model for one dimensional 
heat conduction/mass transport and the identification of a nonlinear dissipation 
mechanism (i.e. damping) in a second order one dimensional wave equation. 
Computational and impleraentational considerations, in particular, with regard to 
supercomputing, are addressed. 
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